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NOMENCLATUEE 


coefficient in velocity expression (eq. (8)), related to upwelling 
index 

argument in confluent -hypergeometric function related to eigenvalues, 
equation (A17) 

coefficient in hypergeometric series, equation (A15) 

mass concentration of biomass 

value of biomass concentration at x = s 

nutrient concentration 

specific heat of water 

upwelled depth 

Coriolis parameter defined by equation (6) 

momentum flux divergence terms equations (2) and (3) 

confluent hypergeometric functions, equation (A16) 

gravitational constant, equation (4) 

energy divergence term, equation (5) 

species flux divergence term, 'equation (7) 

eddy thermal conductivity of water 

defined by equation (A2) 

rate coefficient for production of biomass species i 
argument in equation CA15) defined by equation (A14) 
pressure 

solar radiative flux 
earth’s radius 
time 

temperature 

velocity component in the x direction 



V velocity component in the y direction 

w velocity component in the z direction 

x,y,z Cartesian coordinates and coordinates shown in figure 4 , 
X, Z defined ,by equation (A9) 

a coefficient of temperature in equation (4) 

3 part of Coriolis function equation (6) 

6 defined in figure 3 

e,e' small distances in x 

q defined by equation (16) 

0Q mean latitude of water mass of interest 

X eigenvalues, equation (AlO) 

Uy spectral absorption coefficient 

p water density 

([) defined by equation (A4) 

defined by equation (A4) 

volumetric production rate 
Si earth axial rotational vector component 

Subscripts 
i species 

N nutrients 

o reference condition 

p pressure 

r radiative 

x,y,z partial derivatives with respect to 

V spectral band 


iv 



ANALYSIS OF COASTAL UPWELLING AND THE PRODUCTION OF A BIOMASS 


John T . Howe 
Ames Research Center 


SUMMARY 


An analysis is performed on coastal upwelling phenomena wherein the 
coastal upwelling index derived from weather data is used as input to a set of 
coupled differential equations that describe, in an approximate way, the fluid 
djmamics, radiative transfer in the sea, and the resulting production of a 
biomass. The curl of the wind stress vector is discussed in the context of 
the physical extent of the upwelling structure. An analogy between temperature 
and biomass concentration in the upwelled coastal water is derived and the 
relationship quantified. Theoretical results are discussed in the context of 
observed natural, phenomena, remote satellite or airborne sensing to obtain 
biomass rate production coefficients, and aerospace technology utilization by 
the commercial fishing industry. The intent of the paper is to provide a 
relatively simple framework within which the dominant features of the problem 
are included. The paper is written in a manner which hopefully provokes scien- 
tific interest; but more importantly, is readable and provides useful insight 
to a vast user community of perceptive and extraordinary men and women who 
fish the seas. 


INTRODUCTION 


In the past decade, a great deal of attention has been given to the 
phenomena of coastal -and oceanic upwelling (ref. 1 ), particularly along both 
coasts of North America, and to some extent the coastal water of South America 
and Africa. In the upwelling process, cold, nutrient, rich water is draxm 
from the depths to the ocean surface, where, typically, nitrates (NO 3 ) and 
phosphates (FO^^) and phytoplankton seed are exposed to sunlight in the 
euphotlc zone. In this zone, the phytoplankton grow by uptake of the nutri- 
ents (ref. 2 ), and in turn are grazed upon by the zooplankton and anchoveta in 
a complicated biological process in which the phytoplankton are ultimately 
enriched by ammonium (NH^) excreted by the zooplankton and the anchov.eta. 
Because of the upwelling phenomena, these coastal regions are exceedingly pro- 
ductive in biomass and are particularly rich as a source of protein sought by 
fishing vessels of many nations. Other regions of the oceans are very produc- 
tive fishing grounds as well, such as the shifting boundary of ocean circula- 
tions of the northern and southern hemispheres, locales near the Erben Sea 
mount (2400 km (1500 ml) west of the North American continent), and the 
Emperor Sea mount (about 1600 km (1000 mi) off the coast of Japan) . 

The coastal upwelling phenomena is caused by the Coriolis force that 
arises from the rotation of the earth, the configuration of the landmass, 
including bottom configuration, and is triggered by a wind of sufficient 



strength and duration with a component parallel to the- coast (refs. 3, 4) such 
that the surface water is driven offshore and is replaced by colder water drawn 
from the depths. 

It has recently been reported (ref. 5) that such upwelling phenomena occur 
along the- edge of the Arctic ice pack north of Spitzbergen, initiated by a 
strong wind blowing for a period of about 2 days. Similar phenomena may be 
partly associated with reports of extensive krill production in the Antarctic 
as well. 

Although a great deal of effort has gone into ocean research over the past 
decade, there has been considerable criticism of a disciplinary approach, a 
lack of desired contribution to social goals, and a realization that much of 
marine science lies at some interface between conventional disciplines 
(ref. 6). 

In the context of the present study, many disciplines are involved: fluid 

mechanics, radiative transfer, transport of and finite rate production and 
extinction of species, satellite (and perhaps airborne) technology, and a host 
of others. It is possible to focus on one discipline, such as fluid dynamics, 
to such a detailed extent that the interface with other disciplines cannot be 
described conveniently and in a way that relates to a prediction of the bio- 
mass production. The same may be said of the biological modeling. or any other 
feature relevant to the end object. 

Thus for present purposes, we make every reasonable simplification to 
arrive at some gross behavior of the system and examine what sorts of param- 
eters may be monitored on a global scale that would provide useful information 
to a user community. To the extent that the reader desires, additional details, 
refinements, and features of the problem may be added to the analysis that 
follows. 

It is a pleasure to acknowledge the very helpful comments of Dr. Larry 
Breaker of the U.S. Department of Commerce, NOAA/NESS, and Mr. Robert C. 

Wrigley of NASA-Ames Research Center. 


ANALYSIS 


Figure 1 is an enlargement of a computer-enhanced photograph taken by a 
TIROS weather satellite, ^ which used infrared imagery to display the tempera- 
ture structure of the ocean along the California coast on July 10, 1976. Dark 
tones indicate warm water and light tones show coastal upwelling of relatively 
cold water that originated below the surface. The cold upwelled water pattern 
tends to follow bathymetric features, with fronts extending approximately to 
the edge of the continental shelf (48 km (30 mi) west of San Francisco) and 
closer to shore between Point Arena and Cape Blanco where the bottom 


^Breaker, Larry: Private communication, NOAA/NESS Satellite Field Service 

Station, Redwood City, California, 1976. 
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configuration drops off sharply near shore. A long tongue of upwelled water 
emanates southwesterly from Point Arena, possibly initiated by the underwater 
ridge associated with that point. The tongue is carried southeasterly by the 
sluggish California current. It may also be noted that the waters of the 
San Francisco Bay system are warm (dark toned) , and a warm thermal plume pro- 
ceeds southerly from the Golden Gate Bridge, interrupting the pattern of cold 
upwelled water. 

Such an upwelling pattern is established by a strong wind from the north- 
west, which typically persists periodically for perhaps 3 days at a time during 
the late winter, spring, and summer months. During such a blow, cloud cover 
disappears and enables a satellite photograph to be obtained. As the wind 
diminishes, the cloud cover tends to return, but the upwelling pattern tends 
to persist for a number of days until it is altered by another weather system. 
It should be noted that upwelling "weather" may be quite localized as it is 
influenced by headlands and other features (discussed subsequently) . The 
author has fished off Pidgeon Point in winds and seas that posed an immediate 
health hazard, although a weather update reported light northwesterly winds of 
10 to 20 knots and seas of 0.9 to 1.5 m (3 to 5 ft) from Point Arena to Point 
Pinos. A satellite, of course, can give a detailed panoramic map of the entire 
upwelling system. Usually, one or two such satellite upwelling pictures can 
be obtained every x-7eek, which is sufficient to guide fishing vessels to the 
strongest fronts. As the initial strong winds begin to diminish, commercial 
fishermen begin to venture forth from harbors in a sequence that depends on the 
seaworthiness of their vessels and their individual tolerance to physical 
abuse. 

When the wind stops or changes direction, the upwelling will terminate. 

A typical event occurred on May 6, 1976 in Monterey Bay during a northwesterly 
wind of "small-craft-advisory" proportions. A satellite photograph revealed 
unusually strong upwelling in Monterey Bay, and commercial fish buyers in 
Monterey reported that their facilities were strained to process the very large 
catch of salmon. By May 11, the wind was light from the northeast, a satellite 
photograph revealed no upwelling in Monterey Bay, and at the end of the day, 
Monterey fish buyers reported to the author that no salmon had been marketed 
that day. That situation prevailed until May 14, when a small-craft-advisory 
wind from the northwest revived the upwelling and the fishery. 

There has been an attempt to mark the boundaries of the upwelling regions 
on a mercator projection chart (fig. 2), which is distributed by Humboldt State 
University Marine Advisors (ref. 7} to 11 fishing ports in California. A 
typical chart is shown in figure 2, corresponding to the satellite photp of 
figure 1. 

Some weak boundaries are delineated as well as stronger broad boundaries. 
Whether or not the boundaries are all "fronts" in the sense that the upwelling 
structure is a cell of water that turns downward at the outer boundary is not 
clear. This may be especially true near headlands because of the windshadow 
effect where the headland shelters the region domwind. Upwind of a headland, 
the wind curl (to be discussed subsequently) may be either positive or nega- 
tive, corresponding to stronger or weaker winds offshore. The negative wind 
curl tends to form fronts while the positive wind curl tends to move cold 
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plumes of upwelled water farther offshore (and may produce long tongues of 
upwelled water near headlands). However, downwind of a headland, if the head- 
land provides windshadow, the tendency is toward positive wind curl, and 
boundaries within the shadow may not produce "downwelling" fronts. However, 
if these boundaries extend beyond the windshadow, they may or may not produce 
frontal structure, depending on the sign of the offshore wind curl. 

With respect to figure 1, the weaker "fronts," such as those near Point 
Arena, south, near the San Francisco Bay, and between Monterey Bay and Point 
Conception are "inshore" of stronger fronts. These weak fronts are often 
fished by large numbers of small commercial vessels as a matter of convenience 
and safety — the fronts are roughly 2 hr running distance, or 30 km (~19 -mi), 
from port or shelter. Thus, 4 hr running time, or 60 km (~37 mi), are required 
for a day’s fishing trip. However, productivity of these weak fronts is often 
marginal for the salmon fishery; the weak fronts yield relatively small numbers 
of smaller fish. However, the stronger coastal fronts, which may lie 90 km 
(~60 mi) offshore, tend to be more productive and are accessible to the same 
commercial fleet of small craft, if they know where to find the fronts. It is 
becoming increasingly necessary for them to do so to remain economically 
viable since increased regulation tends to restrict their catch, fuel costs 
rise, and fishing is restricted in both Mexican and Canadian waters. Thus, 
American coastal waters must be fished with greater efficiency for this fleet 
to maintain its utility. The cost is largely due to inconvenience, with 4 to 
6 hr running time, or 60 to 90 km (~37 to 56 mi), each way and quite possibly 
one or more nights spent on the fishing grounds with the catch preserved in 
ice. For the most part, the stronger coastal fronts are accessible to the 
existing commercial fleet of small craft — typically 10 to 15 m (33 to 50 ft) 
in length (a representative vessel and its gear is shown in fig. 3). These 
generally carry modern electronics for communications and navigation — one or 
more two-way radios, radar, Loran, depth recorders and flashers, autopilots, 
and compasses. Moreover, they carry hydraulically operated fishing equipment, 
such as shown at the bottom of figure 3, and have a cruising range 1000 to 
2000 km (621 to 1242 mi). Given the location of the strong fronts periodically, 
these numerous craft can fish them economically. For far-offshore fisheries 
mentioned previously, a fleet of larger, more sophisticated, seaworthy vessels 
is in operation, and can be aided by a knowledge of offshore fishery locations, 
improved weather information, and improved communications from satellite 
technology. 


Mathematical Formulation 

This section is presented primarily for those whose interest is in model- 
ing upwelling structures mathematically. It is a point of reference for the 
analysis that follows. Those who are not familiar with the detailed set of 
hydrodynamic equations may wish to go directly to the next section where cer- 
tain assumptions are made that greatly simplify the analysis, thus sacrificing 
precision in the hope of obtaining a glimpse of the gross features and param- 
eters of the solution. 
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The general equations that describe the water motion (refs. 8 and 9) are 
statements of conservation of mass (eq. (1)), momentum (eqs. (-2), (3), and (4)), 
and energy (eq. (5)):' 
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where the Coriolis parameter or latitude function is 


f (y) = ^sin 0^^ + cos 0^ + gy 


( 6 ) 


and contains the effect of the rotation and curvature of the earth, and F and 
H are the momentum and energy flux divergence term. • Here the usual approach 
of the oceanographer to large-scale phenomena has been adopted that uses 
pseudolaminar models involving constant transport coefficients as free param- 
eters, and the z coordinate is directed upward in a "right-handed" coordinate 
system. In other formations (ref. 10, p. 433), the z coordinate is directed 
doxraward, a convention that will be adopted here, but retains a right-handed 
coordinate system. To these equations we add a species conservation equation 


9c. 9c. 9c. 9c. 

1 , 1 . 1 . 1 
^ + u -r — + V -r— + w 


9t ' “ 9x ‘ ' 9y ' " 9z 
where J is the species .'flux divergence term. 


w. 

J = -^ 
P 


(7) 


Approximate Solutions 

In keeping with the remarks of the last two paragraphs of the introduc- 
tion, we take the local upwelling process to be essentially two dimensional 
and idealize the geometry of the phenomena to the coordinate system shorn in 
figure 4, where the upwelling region is considered to be rectangular for con- 
venience, and the z coordinate is downward. The schematic sho™ in figure 4 
omits the detailed structure of the Ekman velocity spiral for convenience, and 
net transport offshore is simply characterized by an offshore velocity u 
(ref. -11). Moreover, figure 4 shows a distinct frontal interface with the 
ocean. This is generally considered to be associated with negative wind stress 
curl near the front, xjhich is conducive to "frontogenesis," and will be 
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mentioned subsequently. By contrast, positive wind stress curl tends to move 
coastal upwelled water farther offshore where it is mixed with, and diluted by, 
oceanic water. Moreover, for present purposes, we assume that the fluid dynam- 
ics has been established by a wind stress of sufficient duration and the 
upwelling is in the steady state and brings cold nutrient-rich water to the 
surface euphotic zone where it is warmed by sunlight and a plankton bloom 
occurs as the surface water moves offshore. Because these are the phenomena 
we wish to consider, we keep the fluid dynamics very simple. For reference 
purposes, the scales d and 25 are of the order 200 m (660 ft) and 90 km 
(56 mi) (near Point Reyes, for example). 

In figure 4, we might note that at the coordinate origins where x and z 
are zero, and again where the surface waters downwell at the "upwelling front" 
(x = 26), the flow resembles that of a stagnation point (ref. 12). Thus, if 
the alongshore current varies slowly with y locally, the flow is essentially 
two dimensional and the continuity equation (1) is satisfied by 

u = ax , w = -az for 0 < x < 5 (8) 

and 

u = -a(x - 25) , w = az for 6 < x < 26 (9) 

where a can be related to the upwelling index (refs. 3, 4) as follows. The 
upwelling index, I, is the component of Ekman transport directed offshore 
(refs. 3, 4). The Ekman transport is-defined as 

M = Y T X k (10) 

where k is a unit vector directed upward, f is the Coriolis parameter, and 
T is the sea surface shear stress vector resulting from the geostrophic wind 
and is described in references 3 and 4 as 



where_^ is the density of the air, cjj is an empirical drag coefficient, 
and U is the wind vector near the sea surface. The shear stress can also be 
computed from boundary-layer theory of a wind blowing over a liquid surface, 
but the present formulation will suffice. 

Thus, to the order of the depth of the well-mixed photosjmthesis zone 
approximation, the quantity a can be related to the upwellingj index I by 

I = pudj = axpdj^ (12) 

where x should be evaluated at some reference point such as e, and d^ is 
the depth of the well-mixed layer in which photosynthesis takes place and is 
of the order of half a wavelength of the sea or swell. The dimensions of I 
are mass flux of water offshore per unit length of shoreline, in keeping with 
the terminology of references 3 and 4. 
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Thus, the component of the Ekman transport at the coastal boundary is 
related to the geostrophic wind at that boundary. In particular, the V com- 
ponent of the wind parallel to the coastal boundary, directed normal to the 
X, z plane of figure 4 (or down the coast), is of importance. As one moves 
offshore, the variation of that wind component becomes important. Thus, the 
curl of V, or 3V/3x, is an important parameter, expecially if it is diminish- 
ing (negative) with distance offshore. When V becomes zero, one of the 
forces that drives the upwelling vanishes, and the upwelling structure tends 
to be arrested. Thus "a" in equation (8) is related to I at the coastal 
boundary; but the quantity "a" in equation (9) should probably be replaced by 
a quantity a', which is related to the negative curl of V. If the latter 
relationship were well known, the size (26) of the upwelling structure could 
be determined by matching the water velocity u at the juncture of the regions 
described- by equations (8) and (9). The simplest relationship is to say that 
the upwelling structure terminates where V goes to zero. Although inertial 
or momentum considerations may preclude -such simplicity, the relationship may 
be a useful estimate. 

There have .been many occasions when the author left port under very strong 
wind and sea conditions, wondering whether or not the ocean was fishable at 
all, only to find that the offshore environment was much more benign. More- 
over, the front is often marked by a long line, of Portuguese men-o-war, which 
would be blown away by any significant wind. This combination of weather con- 
ditions indicates a negative coastal wind curl, wherein the wind velocity 
becomes negligible at or near the front. In any event, the dimension of the 
upwelled region probably can be related to the wind curl in either an approxi- ■ 
mate or empirical manner. Such a relationship could probably be derived from 
remote sensing data or correlated from surface data. The effect of headlands 
can be determined as well. For present purposes, in a simple analytical frame- 
work in which, to treat the entire problem, equations (8) and (9) are used to 
describe the flow field. 


The energy equation may be written to include the energy flux divergence, 

3T 


pCpH, explicitly as 


3T ^ 9T 3 /. 3T \ 

pucp + pwcp ^ = -JJ (k 


( 13 ) 


where k- is the "eddy conductivity," and q_ is the radiant flux at a depth 


= E 

v=i 




(14) 


Here, q^ is the spectral flux that penetrates the water. It has passed 

through the atmosphere, the air-sea interaction zone (ref. 13), and is the 
unreflected and unscattered residue of the incident flux that actually enters 
the water. A comparison of incident spectral solar radiation with the spectral 
extinction coefficients of seawater (refs. 11, 14) is shown in figure 5. 
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Generally, coastal water absorption is roughly double or triple that of 
oceanic water in the blue, green, yellow, and orange parts of the spectrum. 
Infrared radiation (wavelengths of about 1 pm or greater) is absorbed essen- 
tially on the water surface. 

In equation (14) , the diurnal variation of the spectral solar flux may be 
thought of ‘as "properly averaged" or neglected. In any event, this detail, as 
well as others, can be added to a more precise analysis. For example, multi- 
plying by a sinusoidal time function to simulate the diurnal radiation 

variation would not complicate the analysis excessively. Indeed, such a varia- 
tion could produce results that lend credence to what some commercial fish 
trollers refer to as the "morning bite" and the "afternoon bite." 

In equation (13). we neglect terms involving .9T/3z because energy 
absorbed from radiation is absorbed in surface layers (90% of incoming .solar 
radiation is absorbed in the upper 10 m of the sea (refs. 15, 16)). The 
author has observed a Secchi disc disappear at depths ranging from as little 
as 0.67 to 5 m (‘~2.5 to 17 ft), at the most, in Northern California coastal 
waters. To the. same order (10 m) the surface layers are "well stirred," say, 
to a depth .of ..half a wavelength d;;^ of the local sea or swell. Thus in 

equation (14) is considered to be approximately constant, so that equa- 
tions (13) and (14) become (averaged over the well-mixed layer) 



with Pv^x 1 assumed. Thus,, for present purposes, the entire flux entering 
the water is absorbed.^ For convenience, define 


9 = f (16) 

where e is a convenient constant and could be related to the resolution of 
the remote temperature measurement for practical purposes. Thus, from equa- 
tions (15) and (16) 


dq _ 

q ^ 

2 ^ 


dT 


v=l 


(17) 


^This simple analysis could be embellished by including the Kubelka- 
Munk (ref. 17) equations to account for reflected light, and light that is 
backscattered in depth. 
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which integrates to 


T = 


v=l 


pc adi 


£n n + constant 


(18) 


When ri - 1, we define T = Tq so that equation (18) becomes 

CO 

E V 


T _ ^ ^ " £n q , 1 < q < - 

o pCpad^ - ' - E 


(19) 


Similarly, from equatipns (9) and (13) with the same approximations for 
6/e < q < 2(5/e) 


T - 



pcpadi 


£n 



( 20 ) 


Figure 6 represents the corresponding surface temperature profile (left 
ordinate) for 6/e =500, which is adequate for illustrative purposes. Because 
of the nature of the "stagnation"-type velocity distribution near the coast 
(q = 0) and near the extreme end of the upwelling structure (q ~ 1000) , the 
residence time is large at those extremes and heating of the water is very 
significant compared to that in the intermediate region across the upwelling 
structure. It has been reported^ that as one approaches the upwelling front 
(where the 'streamlines actually turn downward — fig. 4), the temperature may 
rise as much as 9.5 to 17“C. The ocean at the front is also experiencing a 
"stagnation"-type flow as it turns downward and is heated similarly (although 
to a lesser extent because of the lower absorption coefficients) so that a peak 
temperature occurs at the front and tapers off on both sides. This peak tem- 
perature, in effect, marks the offshore front. The asymptotic temperature 
approach to Infinity at the front is a mathematical artifact resulting from the 
approximations made and is not physical. Thus, the solution should terminate 
at some distance, e’, before x = 26, which could be determined from evaluation 
of the dimensionless ordinate and use of the limiting temperature rise cited 
above, or use of' a measured or remotely sensed peak temperature. An alternate 
formulation of the energy equation that does not invoke some of the assumptions 
made, and leads to an analytic solution, is outlined in the appendix. The 
analytic solution is in terms of confluent hypergeometric functions, which may 
provide entertainment for 19th Century mathematicians, but adds little insight 
to the problem. 


The "front" mentioned above is marked in other ways as well. Commercial 
fishermen often look for color change in the water to mark the front, and will 
declare that if they are within 16 km (~10 mi) of such a front, they will find 
it by watching water color (see fig. 5). The reason is that the upwelled 
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water is often green because of the nutrients and presence of a diffuse bio- 
mass, while the ocean water at the front tends to be blue because of a lack of 
these substances. Blue is the "desert" color of the sea. Another mark of the 
front is the flotsam and jetsome, which tend to be carried there by the slug- 
gish current, where their offshore motion is arrested by the oceanic counter- 
current that plunges downward at the front. Another means of identifying the 
front from shipboard is the presence of vast quantities of bait in the water 
near the front. As mentioned previously, the author has observed a line of 
Portuguese men-o-war as far as the eye could see along the front. These small 
creatures that float on the surface are incapable of swimming against the slug- 
gish upwelled current, and in the absence of a strong wind, are simply carried 
to the front by the upwelled water, where they remain. While fishermen search 
for the front, so do the fish. For example, salmon in the euphotic zone are 
thought to have chemoceptors along the length of their body that detect bio- 
gradients. When the gradient peaks, they have found the front and tend to 
disperse throughout the water column along the front. 

Consider now the production of species in the upwelled water. The great- 
est effect of sunlight on upwelled water is concerned with the metabolism of 
living organisms; it provides the energy for the photo synthetic processes of 
microscopic plantlife that all animals depend on for nourishment. The process 
has oxygen as a byproduct, which is of great importance to the respiration of 
animal organisms and in the oxidation of both organic and inorganic substances. 
The life cycle of any system of species is very complicated and in some cases 
can be modeled in a detailed way, such as that mentioned in the introduction. 
That particular model has been programmed by reference 2 in a manner that is 
weakly coupled to the physical world of fluid dynamics and radiative transfer 
as approximate boundary conditions for the biological model. The authors of 
reference 2 acknowledge that the "real world" is not adequately represented. 
Such a detailed biological model could be built into a more complete computer 
code that adequately simulates the fluid dynamics and spectral radiative 
transfer. Such a model would require rate coefficients for each step in the 
biological process. 


It may be noted that the upwelling season lasts several months,. Phyto- 
planckton tend to grow exponentially for, perhaps, 2 days, but then concentra- 
tion levels off because of grazing by the zooplanckton, as they are carried 
offshore by the upwelled water. Thus for present purposes, we illustrate the* 
biological process by simply using the previously mentioned fact that the 
photosynthetic process .of species production ultimately depends on sunlight. 
Accordingly, we assume that the nutrient concentration is essentially invariant 
with x, and is well mixed in the z direction so that the local produc- 

tion rate of biomass per unit volume, is proportional to the nutrient concen- 
tration cjg and divergence of the particular spectral band q.^, (or bands) of 

flux that produces the energy to create a biomass concentration Cj[^. Thus, 
equation ( 7 ) becomes simply (in the quasi-steady state) 


pu 



u. 


1 


( 21 ) 
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where the production rate per unit volume is 


0 ). 

1 


1 N 



( 22 ) 


The biomass concentration Cj^ typifies that which may be inferred from satel- 
lite imagery. Use of equations (8), (16), (21), and (22) in equation (2) 
yields, with the vertical averaging that produced equation (17) ‘ 


apd, dc. 

^ X _ OT]^ 

K . C.^q q 

1 N^V , 

1.. 


which yields, as before 


c. 

X 




apd^ 


in q , 


1 < q < 


(23) 


(24) 


and 


c . 

X 



K - c^^q 
X N ‘v. 

^ 

apd. 


in 




n 2m 


(25) 


which are identical in form to equations (19) and (20) and thus are represented 
by figure 6 (the right-hand ordinate).*^ 


It should be noted that for both ordinates, in accordance with our well- 
mixed assumption, water temperature and biomass concentration are invariant 
with the depth z in the surface layers. The temperature insensitivity to 
depth z for small values of z has been verified experimentally by the 
author, while -at sea, by his thermometer readings. 

With reference to figure 6, the right-hand ordinate for biomass concen- 
tration in the upwelled waters, the biomass increases greatly near the upwelled 
front because of the increased residence time mentioned previously. On the 
oceanic side of the front, there is no biomass — unlike the continuous tempera- 
ture behavior across the front that was caused by solar heating on both sides 
of the front as noted previously. 


Rather, the high concentration of biomass on the upwelled side of the 
front is simply carried doxmward at the front. Three important conclusions 


^An alternate treatment of the species equation, in which additional 
terras are retained, can be formulated like that of the energy equation, and 
would appear much like that outlined in the appendix. 
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may be made. First, the greatest concentration of surface biomass occurs at 
the front in the upwelled water. This is the .place where bait abound and large 
fish feed. It is, of course, the best place to fish. Second, as the upwelled 
front plunges downward, it carries the large biomass, with it so that the entire 
water column along the front also is well supplied. The author has found, by 
commercially fishing the entire water column along the front, that the salmon 
harvest in the upper and middle waters along the front is enormously superior 
in size and number to that away from the front; and that bottom fish along the 
front are also larger and greater in number.- Third, within the approximations 
made, the dimensionless- temperature profile is identical to the dimensionless 
biomass concentration profile, as shown in figure 6. The validity of this 
result could be tested by methods described in reference 15 in which water 
temperature and phytoplankton concentrations are -sensed at wavelengths of 
0.443 pm, referenced to 0.525 ym. Such measurements could yield absolute 
values for Tq and Ci^ as well. The analogy between temperature and biomass 
concentration leads to the relation (eqs. (19), (20), (24), and (25)) 


( 11 %) 
(f - To) 


"p^i‘=N 



(26) 


Thus, the ratio of the biomass concentration increase to the temperature 
increase across the upwelling structure is proportional to the ratio of energy 
absorbed by photosynthesis to the total energy absorbed in the water mass. 


CONCLUDING EEMARKS 


The simplified, analytic mathematical model of an upwelling front that has 
been presented is illustrative of the kind of mathematical study that could be 
done in greater detail by the use of computers. The detailed model could 
include more refined fluid dynamics, radiative transfer, and biological model- 
ing, as well as an improved understanding of the actual wind curl and the 
effect of headlands on the wind curl. It could include a broad geographical 
region (even global) in an efficient way by dividing the region into appro- 
priate sectors and using a parallel processing computer such as the ILLIAC IV. 
The computations could be made on a continuing basis, making use of satellite 
or airborne sensor data to predict biomass distribution and to achieve a better 
understanding of the mechanisms and the kinetic rate' parameters of biomass 
production. The model could reduce the fuel consumption per unit mass of pro- 
tein harvested and increase the protein harvest from the sea. It could dis- 
cover new, untapped fisheries to harvest. It may reveal upwelled regions 
deficient in nutrients, into which nutrients could be air dropped to enhance 
productivity;^ the technology exists to air drop up to 23 m^ (6000 gal) per 
hour. 


^Private communication from C. K. Lombard, Pacific Design, Inc., 
Palo Alto, California, 1979. 
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Moreover, as Ocean Thermal Energy Conversion (OTSC) activities increase 
(especially in areas where conventional fuels are not plentiful) , the use of 
satellites, airborne surveys, and computer simulation of the effects of the 
physical and biological processes associated with an OTEC plant could be deter- 
mined. To a first order, an OTEC plant discharge could be modeled as an axi- 
symmetric upwelling that is not dependent on a wind parallel to the shore, but 
which produces a biomass and a fishery. When the OTEC plant is cleaned of 
biofouling, discharges will affect the biomass. These effects can be deter- 
mined as discussed above, and the rate coefficients for the recovery or normal- 
ization of the biomass can be- Inferred by the same procedure. 

Other uses of multidisciplinary aerospace technology on global problems 
suggest themselves: coastal zone monitoring and management; shoreline erosion, 

monitoring, and protection assessment; the seasonal beach sand distribution 
and budget, including the effect of river dams on shorelines; and the effects 
of artificial reefs and barriers on both the biomass and shoreline protection. 
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APPENDIX A 


ALTEENATE FOBMULATION OF THE ENERGY EQUATION 


From equations (8) , (13) and (14) the energy equation becomes 


3T 3T 3 /, 3T \ 

-5: - -3? ■ 3i- r -aj “Zj ^ 

v=l 


Vv^l 


For constant k, let 


ape. 


and Ko = 


ape. 


so that equation (Al) becomes 


X 


3x 


= K, 


3^T 

3z^ 


+ z 


3T 

3z 


+ K, 


CO 

E 

V=1 


u q 

V V, 




let 


T(x,z) = <j)(x,z) + i|j(z) 

The partial derivatives are then 

^x “ ‘^'x ’ ~ ’ ^zz “ ^zz ’^zz 


^<^x = + ^zz> + + ’^'z) + S 2 

V=1 


y q 

V V, 


“UvZ 


which is now split into two equations 

= ^l^zz + ^^z 


and 


Kl’J'zz + y’J'z = ^2 


V=1 


V q 

V V, 




Let 


4>(x,z) = X(x) • Z(z) 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 
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so that equation (A7) becomes 


vt K,Z" 
xX 1 , zZ' 

X “ Z Z 


(AlO) 


where the primes refer to ordinary derivatives and X~ is a sec of constant 
eigenvalues so that 


xX’ - X^x = 0 

and 


(All) 


KjZ" + zZ' - X^z = 0 
The solution of equation (All) is 


X 



(A12) 


(A13) 


while that of (A12) is 


Z(z) = w(n) (A14) 

where co(n) is a pair of confluent hypergeometric functions, the first of 
which is given by the series 




CO 

p=0 




a(a+ 1) . . . (a + p - 1) 

+ ^ . . . (i + p-i)pi J 


while the second is 


0 ), 


- + ! - I > 2 - I ; ti)' 


(A15) 


(A16) 


and the second series of _iFi is like the first with [1 + a - (1/2)], 
[2 - (1/2)] in place of a, 1/2, respectively, where 


a 



(A17) 


includes the eigenvalues. 

Equation (A8)- can be solved by standard procedures; boundary conditions 
must be worked into the solutions of equations (A8) , (All) , and (A12) in a 
tractable manner, and the eigenvalues must be determined. Evaluating the 
solution of equation (Al) in the manner outlined may prove to be tedious. 
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Figure 1.- Central California coast showing upwelling, July 10, 1976. 
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Figure 2.- Mercator chart of upwelling fronts. 
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TYPICAL COASTAL FISHING CRAFT 



HYDRAULIC FISHING GEAR 



Figure 3.- Typical coastal commercial fishing vessel and its gear. 
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Figure 4.— Physical model simulation. 
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Figure 5.- .Comparison of solar spectrum with optical properties of water- 
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Figure 6 ,- Normalized surface profiles for upwelling zone. 
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